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SUMMARY 


A discussion of the integral relations for flow of the boundary- 
layer type is presented. It is shown that the characteristic laws of 
spread of Jets, wakes, and so forth, can be obtained directly for the 
laminar case and, with the help of dimensional reasoning, for the tur- 
bulent case as well. 

Measurements of the mean velocity, the intensity and scale of the 
turbulent fluctuations, and of the turbulent shear in a two-dimensional 
mixing zone are presented. The results of these measurements are com- 
pared with the mixing-length theories. It is shown that both mixing 
length and exchange coefficient vary across the mixing zone. The theo- 
ries based on the assumption of constant mixing length or exchange 
coefficient are thus in error. 

A discussion of the energy balance of the fluctuating motion is 
given and the triple correlation is estimated. 


INTRODUCTION 


The investigations presented herein constitute one part of a long- 
range research program on the development and nature of turbulent flow. 
The background and scope of this program are outlined in order to show 
why the experimental investigations, of which some simple parts were car- 
ried out some time ago, should be undertaken. 

Theoretical turbulence research can be divided into three main sec- 
tions: laminar stability, phenomenological theories of fully developed 

turbulent flow, and statistical theories. Following Reynolds' work, the 
emphasis of the theoretical research centered mainly on the stability 
problem. The question of whether laminar flow was stable or not stable 
in general became the subject of a great many papers. Both the energy 
method of investigating instability and the method of small perturbations 
were used extensively. It lies in the nature of these investigations 
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that they are analytical theories; that is, the results of the stability 
investigations should lead to results which, without adjusting any con- 
stant, agree with experimental findings. The first marked success of the 
instability theories was G. I. Taylor's discovery of the instability of 
flow between rotating cylinders which was found to be in excellent agree- 
ment with experiment. The theory of stability of two-dimensional flow 
was first successfully carried through for the case of Couette flow. 
suits for the more important cases like plane Poiseuille flow and Blasius 
flow (that is, boundary-layer flow) were doubtful for a long time. The 
investigations of Heisenberg and Tollmien were not considered mathemati- 
cally rigorous, and experimental evidence was not available. Only in the 
last 4 years has the theory of laminar instability been definitely set- 
tled as a result of the work of C. C. Lin (reference l). The results of 
this instability theory for the Blasius flow were found to be in excel- 
lent agreement with experimental results . 

The stability theories deal with small perturbations and linearized 
equations and thus can predict only the lower limit of stability, that is, 
the limit at which a small perturbation will increase in amplitude for 
the first time. The stability theories are unable to predict transition 
and the onset of turbulence. There is little doubt that in both the tran- 
sition region and the fully developed turbulent region the nonlinear terms 
are essential, and this is, of course, the reason why an analytical theory 
of developed turbulence and of transition is nonexistent. Theories o± 
fully developed turbulence are, in general, phenomenological theories. 

Some problems have also been attacked by statistical methods . 

The phenomenological and statistical theories form the second and 
third groups of theoretical turbulence investigations. The phenomenolog- 
ical theory of turbulent shear flow also goes back to Reynolds' introduc- 
tion of the "apparent shear." These theories developed very rapidly, 
following the introduction of the "mixing length" concept by L. Prandtl. 
The main progress was due to the work of Prandtl, Karman, and Taylor. The 
momentum transfer, the vorticity transfer, and the similarity theories 
were developed in rapid succession. The mean-velocity distributions 
computed from the .mixing-length theories were found to be in good agree- 
ment with experimental results at the time. Skin -friction formulas, for 
example, the famous logarithmic law, could be developed on the basis of 
these theories and checked with experiments. Experimental research in 
these years was directed mainly toward finding which of the three main 
theories was most nearly correct. The answer found from the experimental 
results, which consisted mainly of measurements of mean-velocity and mean- 
temperature distributions in boundary layers, channels, jets and wakes, 
was' somewhat disturbing. The flow near a wall was found to agree fairly 
well with results of the momentum transport and similarity theories; flow 
in wake and jets agreed better with the vorticity transport theory. The 
evidence for these various theories has been discussed in reference 2, 
and in papers by Karman (reference 3) and Taylor (reference 4). The 
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shortcomings of these three phenomenological theories have "been pointed 
out clearly "by Karman. In the course of the present report, the authors 
will return to this discussion in some detial. 

It is thought that, at the present time, the mixing- length theories 
have lost much of their, value. As a matter of fact, the main results of 
these theories can he obtained by dimensional reasoning, without the in- 
troduction of a hypothesis of the mechanism of turbulence. The mixing 
length itself, which can now be computed directly from measured quanti- 
ties, is found to be a complicated function of the coordinates. Thus, 
the principal advantage of the mixing-length theories, namely the intro- 
duction of a "simple" length, is gone. It should be mentioned here that 
the realization that dimensional analysis furnishes practically all the 
results of the mixing theories has grown steadily since Karman' s 1937 
paper. C. B. Millikan (reference 5), Mises (reference 6), Reichardt (ref- 
erence 7), and lately Squire (reference 8), have expressed this point of 
view. In the discussion in later parts of the present report this point 
will be taken up again. . 

There remains the third group of theoretical investigations, the 
statistical approach. In the simplest case, that of isotropic ^turbulence, 
first treated by Taylor (reference 9) and subsequently by Karman (refer- 
ence 10), a complete kinematic analysis could be developed. The dynami- 
cal equations of motion could be reduced to a point at which definite 
solutions could be obtained for various specific assumptions. But even 
here there does not exist a physical principle which permits certain se- 
lection from among the various possible solutions. There have been but 
few attempts to extend the statistical theories to more complicated cases 
such as shear motion; Karman discussed these possibilities and gave one 
solution for the case of plane Couette flow. Even more general attempts 
for a statistical theory of. turbulence have been made by Burgers; however, 
no results which can be checked with experiments have been obtained as yet 

Chou has developed a theory of turbulence (reference ll) using the 
equations of correlation starting from Reynolds ' equation for the double 
correlation and proceeding to correlations of higher order. Because of 
the nonlinear term in the equation of motion, it is impossible to obtain 
an equation which includes only correlations of one order. The double 
correlation equation thus involves terms containing triple correlation, 
and so on. The difficulty here is, of course, to find an argument to 
break the subsequent equations off at a certain point. Here again a clear 
physical idea is needed, but as yet has not been found. The application 
of Chou's equation thus involves an arbitrary assumption which, as in pre- 
vious theories, has then to be checked by a comparison with experimental 
results. 

In summarizing, there exists at the present time no satisfactory .the- 
ory of turbulent flow. The concept of a simple mixing length is not 
usable and does not agree with experiment. The recent introduction by Lin 
in an unpublished work on velocity and temperature distributions in 
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turbulent jets, and by Prandtl and Gortler (reference 12) of the assump- 
tion of constant exchange coefficient does not agree with measured values 
of the exchange coefficient, as will be shown in this report. In general, 
the ability of a specific phenomenological theory to predict a mean- 
velocity distribution which agrees well with measurements is not at all a 
proof that the assumptions of the theory are correct. In fact, in the 
example of the plane mixing region which is discussed in detail in this 
report, it will be shown that both the theory based on a "constant mixing 
length" and also the one based on the "constant exchange coefficient" can 
be made to agree with the experimental mean-velocity distribution, in 
spite of the fact that neither the measured exchange coefficient nor the 
measured mixing length is even approximately constant. 

The situation, therefore, appears to be as follows: If only the 

mean-velocity distribution of a turbulent -flow problem is desired, it 
usually can be obtained by dimensional considerations, assuming a reason- 
able curve (e.g., an error function in the case of the mixing region), 
and the determination of one constant from experiment. For an understand- 
ing of turbulent flow in general a study of various mean -velocity dis- 
tributions is of little use. The fields of turbulent fluctuations must be 
studied in detail. At the present time, it is believed that extensive 
experimental research is of primary importance. The turbulent fluctua- 
tions should be investigated in various cases of turbulent flow ranging 
from isotropic turbulence to boundary layers and jets. The main purpose 
of such an investigation would be to determine which characteristics of 
the fluctuating field are universal and which vary essentially from one 
group of flow cases to another (e.g., free turbulence like jets and wakes 
on one hand, boundary layers and channel flow on the other). Eeichardt 
has proposed the development of a theory of turbulence, in an inductive 
way, from measurements (reference 7). However, since the only basis of 
Eeichardt' s approach is mean-velocity distributions, his work is open to 
exactly the same criticism as that applied to a comparison of other theo- 
ries with mean -velocity distributions. Eeichardt' s final equation, there- 
fore, aside from objections which can be raised regarding certain invari- 
ant properties, appears to be merely an empirical Interpolation formula. 

The hot-wire technique is developed to a point at which all fluctuat- 
ing velocity components and all double correlations can be measured. It 
is likely that a knowledge of the triple correlation, that i3, the diffu- 
sion of turbulent energy, is necessary. The development of a hot-wire 
circuit able to measure triple correlations is at present in progress. 

The measurements which are discussed in this report are measurements 
of flow in a two-dimensional mixing region. Earlier measurements along 
the same lines of attack have been carried out on an axially symmetrical 
jet. Measurements in a two-dimensional channel are now in progress as 
are preliminary measurements of a heated jet with essential differences 
in density. 
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To carry out experiments in a field where there is little guidance 
from the theory is somewhat difficult. It is entirely possible that 
after a complete theory is known, it will be found that some trivial meas- 
urements have been made and important ones omitted. This was the case in 
the laminar instability and transition investigations in which many meas- 
urements of "natural" fluctuations in the laminar layer were carried out 
by various observers. After the problem was clearly understood, it be- 
came evident that many of these measurements were of no U 3 e and repre- 
sented merely investigations of complicated effects related to specific 
wind tunnels. 

In order to reduce this possibility as much as possible, it is neces- 
sary to (a) try to define the experimental conditions as clearly as possi- 
ble] for example, make sure that a "two-dimensional" Jet is two-dimensional 
to a sufficient degree, and (b) use as much guidance as possible from 
energy and momentum considerations and from dimensional analysis. 

In a field as complicated as turbulence, even precautions of this 
nature are not always successful. For example, the effect of free— stream 
turbulence upon the isotropic turbulence downstream of a grid has been 
investigated in three institutions; Cambridge, National Bureau of 
Standards, and GALCIT, with different results. The reasons for this dis- 
crepancy are not yet clear. Here, evidently, a quantity which has been 
considered of minor importance and has not been measured cannot be neg- 
lected. This case shows also that a certain duplication of research is 
desirable here rather than superfluous. 

This investigation was conducted at the Guggenheim Aeronautical 
Laboratory, California Institute of Technology, under the sponsorship and 
with the financial assistance of the National Advisory Committee for 
Aeronautics . 

The authors would like to acknowledge the cooperation of Mr. S. 
Corrsin. 


SYMBOLS 


Xi Cartesian coordinates, Xj. = x is the longitudinal coordinate; X— 

axis is defined by the streamline -H- = 0.5; x = 0 corresponds 

to the Jet opening; x 2 = y is the lateral coordinate, increasing 
toward the free stream 
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c a constant determine A by fitting the measured velocity profile to 
the theoretical one 

w magnitude of instantaneous velocity 

W magnitude of mean velocity 

q magnitude of fluctuating velocity 

uq component of w in the xq— direction 

Ui = u component of W in the Xi (^x) direction 

U 2 5 v component of W in the X 2 (=y) direction 

Uq* component of q in the xq— direction (ui* = n*, U 2 1 = v *) 

U 0 free-stream velocity 

p instantaneous static pressure 

p* static-pressure fluctuation 
p density 

p absolute viscosity 

V kinematic viscosity 

6 boundary— layer thickness in general 

A momentum thickness of boundary layer 
7\ microscale of turbulence 
L scale of turbulence 

r shearing stress 

e turbulent exchange coefficient 

l mixing length 

k double correlation coefficient expressing correlation between u* 

and v* at a given point 

Ry double correlation coefficient expressing correlation of u* at 
two points displaced in y-direction 
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R z double correlation coefficient expressing correlation of u' 
points displaced in z-direction 

i mean value of current through hot wire 

i Q current through hot wire when velocity is zero 

e voltage across hot wire 

R mean value of resistance of hot wire 

Ro resistance of hot wire at 0° C 

R a resistance of hot wire at room temperature 

AR R — R a 

S sensitivity of hot wire to flow perpendicular to wire 
Z compensation resistance 
7 theraogalvanometer reading 
G gain of amplifier 
a temperature coefficient of change of resistivity of the hot wire 
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ANALYTICAL CONSIDERATIONS 
Classification of Turbulent— Flow Problems 


The fundamental problems in two-dimensional shear flow, turbulent or 
laminar, can be classified according to the following table: 


Boundary 

Symmetric 

Asymmetric 

Free 

Jet, wake 

Mixing region 

Solid 

Channel 

Boundary layer 



(Couette flow) 


The Couette flow is difficult to realize experimentally. The flow 
between rotating cylinders is similar, but the influence of centrifugal 
forces induces secondary flow even in the turbulent case . 
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The equations of motion for the mean— velocity distribution can be 
reduced to the standard boundary— layer form of the equations for all 
these configurations except the channel flow (and also Conette flow);, for 
which even simpler equations hold. It is therefore useful to discuss 
first the integral re3.ations of boundary flow. 


Integral Relations of Boundary-Layer Flow 1 

To discuss the behavior of the mean— velocity distribution, the inte- 
gral relations for momentum and energy in boundary- layer flow prove very 
useful. The momentum equations, that is, Karmen's integral equations, 
have been used a great deal in considerations of this type. It will be 
seen that the addition of the energy integral equation is very useful for 
a general discussion of the behavior of the mean velocity distribution. 
The boundary— layer equations are: 


pu 


8u 

cbc 


+ pv 


ou _ dp 
dy dx 


+ 


dr 

dy 


(i) 


dpu dpy 
"S + ~"dy 


0 


(2) 


Integrating equation (l) with respect to y between two variable limits 
a(x) and b(x), say, furnishes the momentum integral equation: 


b b 

/ pu 3x ^ + / pv If ^ ~ “ a) Si + T b " T a (3) 

a a 


By transforming the second integral on the left. 


pv ^ dy = 

oy 


pvu 


rV r f r b ^ 

I u j[y - j pvu | + / u dy 

-L dy L J a 4 dx 


n 
I 

a "a 


Hence, equation (3) becomes: 



dy = - [b - a ] ^ 


a P 


dx 



(b) 


i It has recently come to the attention of the authors that a gener- 
alized technique for the application of integral methods for boundary— layer 
calculations has been presented by L. G. Loitsianskii in reference 13. 
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In all cases of free mixing and also for the boundary layer on a flat 


plate = 0, and equation (4) becomes, in the laminar case, 
ox 

b b 

r dpu 2 . r du i 

/ -r — dy = hi r puv I 

J dx I. dy J 

a a 


(5a) 


and, in the turbulent case. 


/ %^dy 

J OX 
a 


— ^pu'v' + puv 


Jb 


(5b) 


Equations (5) will be discussed for specific cases after similar equa- 
tions for the energy have been Witten. 

Multiplying equation (l) with u and integrating furnishes an 
energy integral equation in the form: 

b b b b 

f pu 2 ” dy + / puv ~ dy * - ^ f udy + / u ~ dy (6) 

J dx J dy ox -> J dy 

a a a a 


By partial integration, 

b b 


j' puv dy = r / P v d y = ~ I P u2y 


2 


3y 




b n b 

_ .1 / u * Sled ay 

2 J dy 

a a 


Therefore, 


a 


- b b 

DUV Su ay = i pu a v ,i[,3 ileal 4y 
Oy 2 L J 2 J dx 


Also, 


r b St 

y u s7 dy = 


- -,b b 

I f T ^ U 

< 

3 a 


dy 


dy 


Hence equation (6) becomes: 

b n b 

r d ( p 3 ■. . dp r , 


+ itf-T 


0 2 1 r ln\ 

~ u 2 v | - / t — dy (7; 


J J 

a a 
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If 


op 

< 5 x 


= 0 it follows that, in the laminar case. 


0 /P 3 A , r Zu p , 

N -- V - y ^ L dy 2 


y sv; 


\ 2 


r f du \ 

y “W; y 


(8a) 


and, in the turbulent case. 


8 ( p 3 , 

! — u ; dy = 

dx \2 / 


— r p 2 

jPUU'v' + — U V 


I —7 8u . 

i + j p u v 5? dy 

a a 


(8b) 


The same type of equation can be written for the axially symmetrical 
case and also for the temperature . 


Applications of the Integral Equations 

The use of the integral equations for mixing problems will be illus- 
trated by two examples: the jet and the mixing zone. 

Jet.— The problem is symmetrical; hence the integration from a to 
b is replaced by 0 to b(x). Then 2b is taken as the width of the 
Jet, that is, b denotes the distance from the axis of symmetry where, 
with sufficient accuracy, u = 0. Thus, equation (5) becomes, 

b b 

/ 5 (pu2) dJ = SJ pu2 dy = 0 

o o 

Since 7=0 and u = 0 at y = b, and v = 0 at y = 0 by 
symmetry, equation (8) becomes, similarly, 


d / p 3 , 

s / 2 u dy - 


Hence, 


r _ 8u , 

-I T -r- dy 


J 

o 


oy 


o 


pu 2 dy = M = constant 


(9) 
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for both the laminar and turbulent cases and 


d / p 3 

s/ i u dy ’ 

•o 


r /x \2 

-J >(§) 


1 


1 ° 


p u r v' 


( 10 ) 


g* dy (b) 


In the incompressible, fully developed jet, assume similarity and thus 
put: 

U = U 0 (x) f(q) T) ■ -jf-T 

b(x) 

Equations (9) and (10a) are sufficient to determine u 0 (x), b(x) for 
the laminar Jet. For the turbulent Jet an assumption must be made for 
u’v 1 . By dimensional reasoning, let 


u’v' « u D (x) g(ti) 


Then equation (9) becomes 


pu 0 b J f“(Ti)dq = M 
o 


and equations (10) become 


p du 0 b 
2 *dx 


-U 


f 3 ( "H ) drj = < 


^ flf'Wf dy 


v 


3 f 

puo j g(T)) f 1 (t)) dy 
o 


or, if the constants are omitted for the time being, 


M 

Uo b~- 


where p = constant over the flow field. 


12 


NACATNNo. 1257 


Hence, 


a / 3. n 

— -(u 0 d) ~ 
dx 




< 



in the laminar case 


in the turbulent case 


Laminar 


Turbulent 


u 


o 



u 0 


/v» 


/m y/* 

W 


b 


/v 2 p y2\l/3 

~ I'TV 


b ~ x 


Or, if Uj is defined o.s the initial velocity of the jet, M = pu-;'" 

for a given nozzle size, and the following dimenu ionle as forms can be 
written: 


Laminar 


^ub''l_e_nt 


u i 

b 

x 


/ \i/3 

' l) 


f m \ • 
- ( ..I ) 

\vx/ 


l 


1/3 



= constant 
z 


The constants can be determined if specific velocity and shear distribu— 
tions are assumed. For the velocity distribution a reasonable function 
is, for example. 



X 

Free mi xing z one .- In the symmetrical problem of the 
tegral relations determine the behavior of u Q (x) and of 


jet the two in- 
b (x) . In the 


unsymmetrical mixing region there exists a free stream of constant veloc- 
ity. Similarity downstream means, consequently, that the velocity and 
shear are functions of only one variable rj . Now, however, the mixing 
process is not symmetrical and two equations are needed to determine the 
different rates of spread into the free stream and into the air at rest 
(or into air moving at a different velocity) . Call the width of the 
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region b(l + a) where b is a function of x, and d is a constant. 
The momentum relation, equation (5), furnishes, if v = 0 at y = b, 


or 


Introducing 


D 

/ 

— ab 


dpu 2 


dx 


dy = 0 


— r pu 2 dy - D U 0 -- = 0 
dx- J . 0 dx 

-ab 


= Z 


u = U 0 f (tj ) , 


gives 




Hence, there is obtained the trivial solution b = constant, that is, 
no mixing, or 


f* dT) = 1 


-a 


which determines a. 


The energy equation furnishes: 
b 

1 f d / 3 


il ^ (Dua)dy = -/ T 

OLD — 


-ab 


or, with the turbulent .t. = pU 0 g(T]), since Uo is the only character- 
istic velocity, ; 

1 


* sv ( f f3(i n - i) = -n 


0 2 
Z. TT 

2 u o dx 


/ 


~a 


Uo / 2 

— M- — J f f dT) for lajrni nar flow 
-a 


-pU 0 J gf 1 .. drj for turbulent flew 
—a 


Ik 
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r 


in the laminar case 



constant in the turbulent case 


r 


/ 



v 


in the laminar case 


UqX 


b 


t 


x 


constant in the turbulent case 


The different rate of spread is determj.ned once the velocity profile is 


An extension to the mixing between two streams of different veloci- 
ties is simple. A suitable velocity profile here is 


where <t> (rj ) is the error integral. The same function was used by 
Reichardt (reference 7) and was found to be also the first approximation 
in Gortler's computations which were based upon the assumption of a con- 
stant exchange coefficient (reference 12) . 

Turbulent mixing problems based on the same assumption have been 
discussed recently by Squire (reference 8) . However, it is believed that 
the use of an integral relation for the energy as well as for the momen- 
tum makes the treatment more uniform and allows the determination of the 
power laws for velocity and width of mixing zones in a very direct and 
simple manner. 


A discussion of the energy balance of the fluctuating motion has been 
given by Karman (reference 1$). The relations given here are essentially 
the same. Karman, however, considered parallel shear motion; here it is 
necessary to deal with flow of the boundary— layer type, that is, in the 
present case there is no mean— pressure gradient, but there are two compo- 
nents of mean velocity. 

In order to derive the equations without extensive writing, it is 
convenient to change notation and use Cartesian tensor notation to denote 
the pressure and the components of the velocity: 


known 


f(n) = <Kti) 


The Energy Balance for the Fluctuating Motion 
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Ui = Ui + Ui* p = P + p' 
U 2 = U 2 U 2 * 

U 3 = U 3 t 


Denote the coordinates "by Xj, Xa, X3. Furthermore, the summation conven- 
tion is applied. The Navier-Stokes equation for incompressible motion 
can then be written simply as, 


du 

dt 


1 ^ u i u k 


1 dp 


dXi 


P dx ± 


+ V - 


d 2 Ui 


dXj dXj 


( 11 ) 


Equation (ll) is written for the instantaneous velocity; averages over 
time will be taken only after forming an energy equation. This is the 
essential difference between this equation and the previous energy rela- 
tion given for the mean flow. Multiplying equation (11 ) by u^ furnishes 


1 dujUj 1 dUjUjUk 

2 dt 2 dxk 


1 dpui d 2 ui 

+ u^v 


P dxi 


dxjdxj 


( 12 ) 


Now u-ju-l = Ui 2 + u 2 2 + U3 2 = v 2 . Furthermore, the last term on the 
right side of equation (12) can be rewritten: 


d 2 u-i 1 d 2 uiu-t /dvuN/diu 
v Ui -±- = - v =-=■ - V'( - — - )( - — - 


dxjdxj 2 dxjdxj v - dx 1 ^ dx j 

Hence equation (12) becomes 


1 dw£ 1 dw 2 u 

2 dt 2 dxk 


Now 


w 2 = U1 2 + U 2 2 + 2(UjUi f + U 2 u 2 ’) + Ui ,cl + u 2 * + u 3 * 2 

= W 2 + 2(UiU 1 * + U 2 u 2 ' ) + q 2 


(13) 


1 dpui 1 d 2 w 2 

p dxi 2 dxjdxj \ dxj> \dxjZ 


Take the time average of equation (13) > and use the following facts: 
(a) The mean motion is assumed steady. 
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(b) The terms r,Cirt^ulai.’lg only mean -speed terms cancel each other. 

(c) The mean value of any fluctuation component vanishes. 

There is then obtained 


i 7 dUx ~ 7 <^J 2 ^ ^ u k<I 2 1 cip’up’ 1 d 2 q^ 

Ux*U k ’ + u 2 »u k 3 - -r^- + - v t— ~— 

6 x k dxk 2 dx k p oxi 2 dxjoxj 

where 


v *£ f-Hi ( 1 i 0 
dxj <3xj 


3u k q 2 _ dU x q 2 + dU 2 q 2 + 


c)x k 


dxx 


6 x s 


ox k 


Fquation (l4) can be simplified somewhat by using the boundary— layer 
approximations. It must be kept in mind, however, that the boundary-layer 
conditions apply to the mean speed and the derivatives of mean quantities 
only, that is, 


Ug, 

Ux 


0 



if 6 denotes a boundary— layer thickness, but 


also. 


however, 


dq 2 

ox 2 


du., * 
5x 2 - 


Ux‘ 


0(1) 



O', 


'St* 
V A / 


i - 


oq 2 

OXx 



\ 

) 


where A denotes the so-called microscale of turbulence and is discussed 
in detail later. By I’.sing this reasoning, equation (l4) becomes: 


ux'u 2 f 



1 du 2 f q 2 

2 ' 6 x 2 


i dpiU2t + 1 v °fi! _ v ^±L ( 15 ) 

P dx 2 2 ox 2 2 Sxj Sxj 


The order of magnitude of the terms of equation ( 15 ) can new be compared 
in a way similar to that used in Karman : s paper: the flew represented by 
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equation ( 15 ) contains one characteristic velocity U and two character- 
istic lengths 6 and T*. The fluctuating quantities will, as before, be 
referred to U. Thus, rewrite: q 2 ~ U 2 , and so forth, and the follow- 

ing expression is obtained: 

3 3 3 2 23 

U U U v U „ IT 

ai — + a 2 — = - a 3 — + a 4 v _ - a 5 v -5 

6 6 6 6 2 * 


where the have the character of correlation coefficients and depend 

on g = eay. In comparing the order of magnitude of the terms in 
equation ( 15 ) the following distinction must be made: 

(a) Compare the order at a given q, at different downstream posi- 
tions . 


(b) Compare the order at different values of T| across the mixing 
region. 

The discussion for item (a) is equivalent to Karman's discussion. 
First, 

' u 2 u 2 

v a 4 *2 « va 5 ^2- 

since S » A, and it is very unlikely that a<i is much different from 
a 5 . The rest of the terms are of the same order if the fundamental rela- 
tion between 7k and 6 is satisfied (reference 14), that is, if 



where a depends on the coefficients aj_ and thus aU ~ q. That this 

type of relation is satisfied will be shown later from experimental re- 
sults in a mixing zone . 

To compare the order of magnitude across a mixing zone or boundary 
layer, as indicated in item (b), requires the knowledge of the mechanism 
of turbulence or experimental measurements . This corresponds, of course, 
to the fact that the mean— energy and momentum equations cannot furnish 
the velocity profile but lead to general relations for the x(downstream) 
variation of the characteristic quantities after a velocity distribution 
is known. The general relation is between A and 6 . 
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From measurements , as reported later, the following picture is ob- 
tained: 


—7- 7 dU x 
Ul U2 5 ^ 


measured ("Production" term) 


3 ^ U2 '^ 


1 _a_ 

p 8x 2 


p»u 2 * 



unknown ("Diffusion" term involving 
triple correlation) 

unknown (same general nature) 


measured 




known approximately (since A(rj) has 
been measured) 


It is evident from the present considerations that measurement of the 
triple correlation is extremely desirable. In the following sections of 
the report a discussion of the measurements for the case of the plane mix- 
ing zone is presented. The results of these measurements will then be 
treated in the light of these considerations. 


EQUIPMENT AND PROCEDURE 
Wind Tunnel 


The investigation was carried out in the wind tunnel shown in figure 
1 . The tunnel is especially designed for investigations of two-dimensional 
nature because of an 8:1 aspect ratio (60 by 7.5 in-). The turbulence 
level is controlled by a honeycomb and seamless precision screens, followed 
by a 10:1 contraction. The screens have 18 meshes per inch and a wire di- 
ameter of 0.018 inch. The honeycomb consists of paper mailing tubes, 6 
inches long and 1 inch in diameter. 

The jet emerging from the contraction (fig. l) is allowed to mix ou 
one of its boundaries with the still air. The other boundaries are solid 
walls. It is especially important that the boundary opposite to the mix- 
ing zone under investigation be closed. This is achieved by means of a 
plate-glass surface which reduces influences of any draft in the room on 
the half jet and improves the two-dimensional character of the jet. 
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The tunnel is operated by a 62-horsepower stationary natural- gas 
engine, normally operating at a fraction of its rating, which drives two 
eight -blade fens. The speed is controlled remotely by means of a small 
electric motor which drives the throttle through a gear and lead-screw 
system. The velocity range of the tunnel is about 5 1° 40 meters per 
second. The experiments were carried out at a speed of 18 meters per sec- 
ond. 


Traversing Mechanism 

Figures 2, 3 } and 4 show the types of mechanism used during the ex- 
periments. Traversing in the longitudinal direction (x -direct ion), par- 
allel to the mean flow, was done by a hand-operated screw (fig. 2). In 
the lateral direction (y-direction) the hot-wire carriage was moved by a 
small 6— volt direct— current motor (fig. 3). A revolution counter indi- 
cated the position of the hot wire within ±0.02 centimeter. 

The hot-wire carriage (fig. 3) consisted of a Z-shape brass bar which 
could be rotated in a horizontal plane by means of a small gear and lead- 
screw system. Means were provided for reading the rotation angle. The 
hot-wire holder consisted of ceramic tubing fitted in a short brass tube 
mounted on the top of the Z bar. The mounting was so constructed that 
the hot-wire holder could be rotated in horizontal and vertical planes. 
This arrangement was necessary for measurements of v* and correlation 
coefficient k. 

For the measurements of the % correlation, the traversing mechan- 
ism shown in figure 4 was used. One of the hot-wire hol.ders was station- 
ary, the other was movable by means of a hand-operated screw. Care was 
taken to maintain the two hot wires parallel. The initial distance be- 
tween the wires was measured with an ocular micrometer. The accuracy of 
the mechanism was ±0.006 centimeter. 


Hot-Wire Equipment 

The hot-wire apparatus used during the measurements is described in 
detail in reference 5* 

Mean-speed measurements .— A half-mil (0.0005 in.) platinum wire of 
1— centimeter length was used. The measurements were made by the constants 
resistance method. This method has the advantage of keeping the wire 
sensitivity constant throughout the velocity profile. 

It must be kept in mind, of course, that the hot wire measures the 
absolute magnitude of the velocity. However, as will be seen later, it 
is justifiable to assume that the lateral component of velocity is very 
much sma l ler than the longitudinal component throughout the main portion 
of the half jet. 
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The transverse component of the velocity was also measured at 
x = 5^.3 centimeters* The method was as follows : Two identical hot wires 

were placed at a distance of 15 centimeters, one "behind the other in the 
downstream direction. The velocity difference Au was measured across 
the mixing zone. Since 

» 

Au du _ dv 
Ax ~ ox dy 

the graphical integration of the function 


gives the distribution of the transverse component of the velocity for a 
constant value of x. 

Turbulence me a surement . — For the investigation of turbulent fluctu- 
ations, 0.00024— inch Wollaston wire was used. The wire was soft- soldered 
to the tips of fine sewing needles after the silver coating had teen 
etched off. The longitudinal component of the fluctuations was measured 
with a single hot wire, the lateral component with a bi-plans X-type 
meter composed of two hot wires . The measurements were carried out by 
conventional methods . 

Measurement of the double correlation coefficient u * - v --.— The method 

u 1 v* 

used is similar to the one proposed by H. K. Skramstad (reference 16) . 

The correlation coefficient was measured with the same X— type meter as 
was used for v' measurements. The method consists in putting the out- 
puts of the two wires separately through the amplifier, and also the sum 
(or difference) of the outputs. It is essential to make a correction for 
the difference in sensitivity of the two wires. The method of correction 
is given in appendix A. This method was worked out and used by P. Johnson 
at the California Institute of Technology during similar measurements on 
a flat plate . 

An alternative method of measuring the double correlation coefficient 
was also used; this involved photographic recordings of correlation fig- 
ures from the screen of an oscilloscope. This type of measurement was 
first made by Reichardt (reference 7). The procedure is as follows: u' - 

and v' -meters are mounted next to each other as close as possible. 

The signals from the meters are amplified by two separate amplifiers 
which must, of course, have the same frequency and phase-shift character- 
istics. The separate outputs of the two amplifiers are then connected, 
one to the horizontal and the other to the vertical deflection plate of 
the oscilloscope. The gain of the oscilloscope amplifiers is set in such 
a way that the amplitudes of the horizontal and vertical fluctuations 
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become the same. The crosses on the pictures (figs. 5 to 11) are obtained 
by setting the gain of first the horizontal and then the vertical oscillo- 
scope amplifier to zero; thus a check for the amplitudes of the fluctua- 
tions is obtained. The figures were photographed with an exposure time of 
15 seconds. The correlation coefficient is calculated from the equation 

k = - a 2 - b 2 
a 2 + b 2 

\ 

where a and b are the major and minor axes, respectively, of the cor- 
relation ellipse . 

Measurement of micro scale turbulenc e.— This involved measurements of 
correlation between longitudinal velocity fluctuations at two different 
points. Two parallel 0 . C0024—inch wires were used, one stationary and 
one movable. From the outputs of the wires, R y was computed. The micro- 
scale of turbulence then was obtained by fitting a parabola to the upper 
part of the R y correlation curve. 

Correlation figures were also obtained on the screen of an oscillo- 
scope. 

RESULTS AND DISCUSSION 
Velocity Distributions 


Careful velocity measurements and visual observation of the velocity 
fluctuations on the oscilloscope show the following facts: 

(a) The boundary layer at the mouth of the jet is only 0.1 centimeter 

thick and is laminar. 

(b) The free boundary layer is laminar from x = 0 to x = 6 centi- 

meters. The transition from laminar to turbulent flow at 
x = 6 centimeters can be seen clearly on the oscilloscope. 

(c) Fully developed turbulent— velocity profiles are obtained only for 

x > 30 centimeters, approximately. 

Figure 12 shows the mean— velocity distributions across the mixing 
zone at different stations downstream: at x = 10, 30, 54.3, 75, and 90 

centimeters. For convenience in the figure the x~axis is chosen parallel 
to the free-stream direction. (This is Tollmien's notation; in all other 

figures the x-axis corresponds to the — = 0.5 streamline.) The 
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boundaries of the mixing region were chosen arbitrarily, defined by the 


that for x > 30 centimeters, where the fully developed turbulent veloc- 
ity profiles exist, the spread is linear. 

In order to demonstrate more explicitly the difference between the 
velocity distributions in the partially developed and the fully developed 
turbulent mixing regions, figure 13 shows the measured velocity profiles 
plotted on a nondimens ional scale. It is immediately obvious that while 
the velocity profiles measured at x = 90; 75; 65 , and 5^*3 centimeters 
(the fully developed region) follow the same curve, that is, they are 
similar; at x = 20 the measured profile deviates appreciably from the 
fully developed type . 

The results of the lateral— velocity-component measurements are given 
in figure 14 . The measurements were carried out only up to | = -0.6, 
because at larger values of | the front hot wire was in the region 
where u is no longer large compared with v; therefore, the assumption 
that the absolute magnitude of the velocity indicated by the hot wire may 
be considered equal to the velocity component, u, no longer holds. 

By use of Gortler's theoretical values of the lateral mean-velocity 
component, which are convenient and accurate enough for computing a co:n- 
rection of this type, u/U 0 was calculated from the measured absolute 

magnitude of the velocity. It was found that, in the region 00 < | < — 1.4 
the error involved in neglecting this correction is of the same order of 
magnitude as the experimental scatter. This justifies neglecting the dif- 
ference between the absolute magnitude of the velocity and its longitudinal 
component for 00 < | < —1.4, and the hot-wire measurements may be con- 
sidered to be correct in this region. 

One mean— velocity distribution, that at x = 5^.3 centimeters, was 
measured with a total— head tube in order to see how closely the results 
agreed with the hot-wire measurements . The results are plotted in figure 
17; the agreement seems to be satisfactory. The effect of large velocity 
fluctuations on the high total— head— tube reading at low velocities can be 
seen clearly in the figure . 

On comparing the results with measurements given in the papers of 
Tollmien (reference 18) and Cordes, the agreement seems to be satisfactory; 
Beichardt’s results, however, diverge appreciably. The value of 0 


ments (reference 19) is 11.99; Cordes gives 0 = 11.95 (reference 15) as 
compared with 12.0 obtained with the present measurements (fig. 14), in 



= 0 . 10 . 

u o 


It can be seen 



\ 

in Tollmien’ s notal obtained by the Gottingen measure— 
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each case the value of a being chosen for the best fit to Tollmen’ s 
curve. Fitting Gortler ’s theoretical curve the present measurements give 
a = 11.0 (fig. 15 ); while Reicharat’s value of a is 13-5 (reference 12 ). 
Reichardt noticed the discrepancy between his and earlier measurements 
and explained it by the fact that early measurements did not go fax* enough 
away from the jet opening. It is believed, however, that three-dimensional 
effects may have disturbed the flow conditions during Reichardt’ s measure- 
ments. The two -dimensionality of the jet used during the present measure- 
ments was checked and confirmed . 

In a comparison of the measured velocity profiles with the existing 
theories it was found, as mentioned in the Introduction, that by an ap- 
propriate choice of the constant 0 both the theory based on constant 
mixing length and the one based on constant exchange coefficient can be 
made to agree with the velocity distribution obtained by measurement. 

Figure 14 shows that for 0 = 12.0 the agreement with Tollmien’s veloc- 
ity profile is fairly good; while for 0 = 11.0 Gortler’s velocity pro- 
file is better approximated (fig. 15 ) • As a matter of fact, an error 
integral curve (Gortler 1 s first approximation) gives a reasonable agree- 
ment with the measured values as shown in figure 15 * 

In order to avoid this arbitrariness in the choice of 0, a conven- 
ient parameter was chosen that could be completely defined from the meas- 
ured velocity distributions at every value of x. Such a parameter is #, 
known as the momentum thickness in the boundary -layer theory, and is de- 
fined at a given value of x as 


Therefore, in presenting the measured data in the figures subsequent to 
figure 15 , all quantities are plotted against y/d. 


Measurements of the longitudinal and lateral components of the veloc- 
ity fluctuations were carried out at three stations: x = 30 centimeters, 

x = 54.3 centimeters, and x = 75 centimeters. Figures. 16, 17, and 18 
show the distribution of the relative velocity fluctuations or turbulence 
levels, u‘/ u and v ? /u, at the three stations. It can be seen from 
these figures that 

(a) The lateral distribution of u’/u remains very nearly the same 
for different distances downstream of the jet opening; the 
same is true of v*/u* 



Turbulence Level 
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(b ) Both turbulence levels reach a maximum at the outer edge of the 
mixing, zone. 

It is interesting to note that Corrsin’s measurements in an axially 
symmetrical jet (reference 20) show similar distributions of both u*/u 
and v* /u. Of course, the turbulence levels in the high-velocity region 
of each cross section are much higher because a fully developed free- 
turbulent jet has no free— stream flow; that is, it has no adjacent stream 

with — - - = 0 and with a turbulence level which is small and due to the 
<iy 

tunnel only and hence, independent of the mixing mechanism. Nevertheless, 
the general character of the distributions is similar in the two cases. 

In order to get a better picture of the velocity fluctuations them- 
selves, u* and v’, their distribution is plotted in a dimensionless 
form, u*/Uo and v ! /Uo, in figure 3.9* It is seen that the va3.ue of 
u 8 /Uo is considerably higher across the mixing zone than the value of 
v ! /U 0 . Furthermore, the v* fluctuations seem to reach maximum value 
somewhat closer to the free stream than do the u* fluctuations. Also 
it should be noted that the fluctuations are somewhat lower for x = 30 
centimeters, because the region here has not yet reached the fully devel- 
oped turbulent state. 

It should be mentioned that no effort was made to correct the hot- 
wire measurements in regions of high turbulence level. The nature of 
these corrections is discussed briefly in reference 20. 


Correlation Coefficient Measurements 
Figures 1 6, 17. and 38 show the distribution of the double correla— 

U. * V " \ 

tion coefficient k = -7- -- at the stations x = 30 centimeters, x = 54-3 

centimeters, and x = 75 centimeters. A comparison of the three distri- 
butions is given in figure 19 where the solid line indicates the average 
values of the correlation coefficient at the three stations. 

The correlation coefficients on the free— stream side of the mixing 
zone seem to be somewhat high. This probably is due to the fact that the 
measurement of the wire sensitivities Sj. and S 2 is more difficult 
here because of the one-sided fluctuations, which will be mentioned later. 
Furthermore, the influence of the ratio S 1 /S 2 on the value of k is 
more pronounced in this region (appendix A) . 

As mentioned earlier, an alternative method of measuring the correla- 
tion coefficient was also employed. These measurements, however, are 
less accurate. This becomes particularly evident near the free stream 
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(fi G s. 6 and 7) ■ It can be seen that calculation of the ^^eUi rse"T 

ficlent (that U, by measuring the major and minor «es of the 

from figures 6 and 7 is urreliable From f 'gures^. " M oumparei wlt h 

correlation coefficient ^ tbe lOTble hotwire method. 

to toother hand, the pictures show the interesting way in which the 
correlation between »• and V changes acrose the 
5 is taken in the free stream where k = 0. It is worth noticing _ 
the circular spot is formed by individual arcs of more or less ciiou^ar 
shate. This can he seen even better in figure 20, where on the right 
si<?e a picture of short exposure shows the shape of such an arc. T 
ttof although h - 0 in the free stream, 
free stream are not entirely of random character, vhich xS c. well Un- 
characteristic of low-turbulence tunnels. The picture inf igure 1 
taken in the isotropic turbulent field behind a screen. Again *-v, 
but the picture on the right shows that here u' and v ar. trul„ 
random character. 

Figure 6 was taken at the free-stream edge of the mixing zone. It 
is seen that the circular shape of the correlation figure is someuhat 
distorted, an indication that the correlation coefficient is 
zero. It is to be noted that this distortion occurs in only one <3.irec 
tion In figures 7 and 8 the distortion is even more pronounced. 1-gUx 
9 shows the highest correlation; while fig-ores 10 and 11 were taken on 
the outer edge of the mixing zone (that is, the zero— velocity edg ) 
the correlation coefficient decreases. 

Corrsin (reference 20 ) also measured the correlation coefficient at 
one section in the axially symmetrical jet. The maximum value was -0.42. 
The explanation for the difference between that maximum value and the one 
presented here is not evident. 


Turbulent Shearing Stress 

From the measured correlation coefficient and the components of the 
turbulent fluctuations, the turbulent shear can be calculated easily. 

Figure 22 shows a comparison of the measured shear and that obtained 
from Tollmien's and Gortler’ s theories. The method of computing the 
shear distributions from the velocity profile of Gortler was given Dy 
C. C. Lin (appendix B) . The same method was used for calculating the 
shear from Tollmien's velocity profile; of course, this can be carried out 
analytically since the velocity distribution i3 given in a closed form. 

It is seen that the measured shear is considerably smaller than the 
theoretical values of Tollmien and Gortler. The difference between the 
maximum values is about 25 percent . The difference on the outer edge of 
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the mixing zone is even larger and is, of course, due to the fact that 
both theories assumed fully turbulent flow aci-oss the entire width of the 
mixing zone. This assumption was found to be incorrect. (See oscillo- 
grams in reference 20 . ) 

The calculation of the shear distribution from the measured velocity- 
profile was also carried out. In evaluating the constant of integration 

(appendix B) the condition that t be maximum where — — = 0 was used 

dy 2 

instead of the boundary condition at the outer edge of the mixing region; 
that is, t = 0 when = 0 . \ 

The reason fox- choosing this was that at the outer edge of the mix- 
ing region the simplified form of the flow equation used in the analysis 
does not hold; it was thought, therefore, that the evaluation of the con- 
stant of integration from conditions at the outer edge is not justified. 

'2 

The condition that t has a maximum at = 0 appears physically 

dy2 

sound. (The mixing— length theories lead also, of course, to this condi- 
tion.) It should be pointed out that by this method of calculation the 
shear has poeitive values at the outer edge of the mixing zone. 

Figure 22 shows a good agreement between the calculated shear from 
the measured velocity profile and the measured shear. 

From the measured shear and the velocity gradient, the mixing length 
and the exchange coefficient distributions across the mixing zone were 
computed (fig. 23). The irregular shape of the curves is due to the dif- 
ficulty of obtaining the slope of the velocity profile graphically. 

It is seen from figure 23 that neither the assumption of constant 
mixing length nor that of constant exchange coefficient holds throughout 
the mixing region. 


Microscale of Turbulence 

t 

Measurements of correlation "between the longitudinal fluctuations at 
two different points at various positions in the mixing region are pre- 
sented in figures 24 to 29. The measured values of R y were plotted on 

a rather large scale, so that the parabola defining the microscale of 
turbulence could be obtained with reasonable accuracy. Figure 24 
shows the Ry distribution at different lateral positions for x = 30 
centimeters. It is seen from this figure that the parabola corresponding 
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to A = 0.24 centimeter is fairly veil defined by the measured points 
and that the parabolas corresponding to A = 0.21 centimeter and 
~h = 0.27 centimeter fall above and below the measured points, respec- 
tively. Thus the determination of ^ by fitting the By curve in the 
neighborhood of the origin can be obtained within about lO percent ac- 
curacy . 

Figure 24 is reproduced in figure 25 in order to show more explic- 
itly the interesting result that the microscale of turbulence is con- 
stant across the mixing region. This result is confirmed further by 
measurements at the cross sections x = 5^*3 centimeters and x = 75 
centimeters (figs. 27 and 28). 

Considerable error in measurements was introduced at points close 
to the outer edge of the mixing zone. At these points, the large lat- 
eral component of the velocity caused the outer wire to fall into the 
wake of the inner one as the two wires were placed close to each other. 


Figure 3 Q shows this effect. It is seen that, while for |- = - 0.97 the 
effect shows up only when the wires are 0.025 centimeter apart, for 

X = —3. h8, the error becomes appreciable at a wire distance of 0.10 

centimeter. No effort was made to correct such error; measurements at 
such points were disregarded since any correction at these large turbu- 
lent levels is doubtful. 


As a matter of interest, some correlation pictures were taken at 
the point x = 75 centimeters, 7 = 1.99 f°r wire distances d = 0.027, 

V 

and d = 0,192 centimeter. (See figs. 31 and 32 .) The pictures on the 
left-hand side were taken with a 15— centimeter exposure, the ones on the 
right were instantaneous pictures. The calculated correlation coeffi- 
cients from these pictures are Ry = O.985 and Ry = 0.89 as compared 
with O.99O and 0.902 obtained by direct measurements. 

The correlation coefficient R z , expressing the correlation be- 
tween the u* fluctuations at two points displaced in the z-direction, 
was also measured at one point of the mixing zone to check the assump- 
tion of two-dimensional flow. Figure 33 shows that the By and Rz 
distributions -are identical within the range of scatter. 


Scale of Turbulence 

From the measured Ry distributions, the magnitude of the scale of 
turbulence could be estimated. This wa3 done by fitting an exponential 
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function to the Ry distribution, as shown in figure 34. It was found that 
for a given value of x the scale decreases somewhat from the free- 
stream side of the mixing region toward the outer edge. The longitudinal 
distribution of the scale at approximately the same values of y/d shows 
a linear increase with x. (See fig. 35-) This figure also shows the 
variation of the microscale of turbulence with x. This variation can 
be approximated by a parabola (straight-line approximation would give a 
finite value of x = 0) . 


Energy Balance Estimated from the Measured Quantities 


The energy equation in a two-dimensional mixing region has the form 
(from equation (15) under Analytical Considerations): 


I du 
p dx 2 





8 g q 2 

dx 2 2 


duj 1 8uiJ_ 
dxj “ 8xj 


The first term represents the production of turbulent energy and can be 
calculated from the measured shear stress and velocity profile. The 
third term can also be calculated from the measured fluctuating veloci- 
ties and can be shown to be negligibly small compared with the other 
quantities. The fourth tern represents the energy dissipation due to 
viscosity. The magnitude of this quantity can be approximated from the 
measured microscale of turbulence. Each term of the dissipation 


D = v 5V 1 
oxj 8xj 


can be expressed in terms of a microscale. If the turbulence is iso- 
tropic, the dissipation can be written in terms of a single microscale 
A: 

,2 

D = 15 V (reference 9) 

A 


A being a function of Ry. During the present measurements, only this 
^ was measured, and although the turbulence in the mixing region is not 
isotropic, as a first approximation the dissipation was calculated from 
the foregoing relation. 

Figure 36 shows the turbulent energy production and dissipation as 
a function of y / The difference of the two function gives, according 
to the energy equation, the diffusion term. It is interesting to notice 
that all three energy terms reach their maximum value at about the same 
region, namely, close to the inflection point of the velocity distribu- 
tion. 
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From the diffusion term, the distribution of the energy transport 
through a unit surface can he estimated. It is seen from figure 37 that, 
in the middle of the mixing region, this term is small compared with the 
double correlation coefficient and becomes gradually larger on approach- 
ing the two edges of the mixing region. At the edges, however, no con- 
clusive remarks can be made as to the relative magnitude of the triple 
correlation, since all simplified assumptions made in the evaluation of 
the diffusion term fail there. 


CONCLUSIONS 


Measurements of the fie.ld of fluctuating velocities in a turbulent 
mixing zone show that both the mixing length and exchange coefficient 
vary across the mixing region. The theories of Toll mi en and Prandtl— 
Gortler which assume constant mixing length and constant exchange coeffi- 
cient, respectively, are thus based on invalid assumptions. It is typical 
of phenomenological theories of turbulence that the mean-speed distribu- 
tions derived from such theories can be brought into fair agreement 
with experimental measurements. This is found to be the case here also if 
the results of Tollmien and Prandtl-Gortier velocity distributions are 
compared with the measured distributions. Conclusions regarding the phys- 
ical significance of theories of turbulence must be based on comparison 
with experimental investigations of the field of fluctuating velocities. 

The discussion of the energy and momentum integral relations for the 
mean and fluctuating motion shows that the over-all characteristics of a 
turbulent mixing process can be obtained by dimensional reasoning without 
any assumptions for the physical mechanism of turbulent motion. 

The following results of the measurements of turbulent field appear 
of general importance: 

(a) The microscale of turbulence A was found to be constant across 
the larger part of the mixing zone. 

(b) Karman's fundamental relation between ~K and the scale of tur- 
bulence L was found to hold. 

(c) The double correlation coefficient reaches a maximum value of 
approximately -O .55 near the inflection point of the mean-velocity 
distribution. 

(d) The three energy terms — production, diffusion, and dissipation 
of energy — are found to have a maximum value in the middle of the mixing 
region . 


Guggenheim Aeronautical Laboratory, 

California Institute of Technology 
Pasadena, Calif., July 24, 1946 
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APPENDIX A 

METHOD OF CORRECTION FOE THE DIFFEEENCE IN SENSITIVITY OF TWO HOT WIEES 
USED FOE C OEEELAT ION--C OEPFIC PENT MEASUREMENTS 


The equilibrium equation of the hot wire is given in the form (ref- 
erence 21) 

Sifm = A 1 + B* */u 
AR 

where A* and B 1 are constants . 


With 


A = — - and B = — - 


oEr 


oE 0 


the equation becomes 


.2. 


— = A + B /u 
AR 


CD 


With i constant the differential of equation (l) is 


i 2 Ra 

(AR) S 


dR = 


du 


- B Su 
2 u 


since i~ 2 = ^-^5 by definition 
C R 


= i (if£ _ a) 5 ? 

2 \AR J u 


1 IL f j_2 _ i 2^j 

2 AR 


VT 


or 


du 


u 


2B al 2 dE 
ARE(i 2 - i 0 2 ) 


Let e = id E and S(iiR) = wire sensitivity = — i2_l 


IB, 


then 


du _ 2e 
u S 


( 2 ) 
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On impressing the voltage e across the amplifier input the thermogal— 
vanometer reading will be 

e 2 = G 2 — = D 2 ? (3) 


where 



If U is considered as the component of the mean velocity perpendicular 
to the hot wire, and the effect of the flow parallel to the wire is neg- 
lected, the equations for uniform, unsteady flow across two wires, as 
shown, are found in the following manner: 


Ux = u sin a 

dUi = sin a du + u cos a da 
Let du = u*, and u da = v* 



dUx = u* sin a + v* cos a 

dUi u* v* ' 

— - = — + — cot a 

Ux u u 

U 2 = u sin (0 — a) 

dU 2 = sin (0 — a) du — u cos (0 — a)da 
= u‘ sin (0 — a) — v* cos (0 — a) 


w 


= — - - V - cot (3 - a) 
U 2 u u 


or for the case 0 = 2a, 


dU 2 u’ v’ ^ 

U 2 u u 


( 5 ) 


For the case a = 45°, equations (4) and ( 5 ) become 
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dUi _ u ’ v f 

U x u u 

dU 2 u» v* 

U 2 u u 


From equations (2) and (6a), 


From equations (2) and (bb) 


I (Hi) 4 2 aLiL + (il 

4 LVu.y u2 v U 
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5 , ft fAff _ 2 + (Zl 


° 2 " 4 LV U V " 

From equations (2) and (6a) and (6b) 


N2 1 

u/ ! 


M®i + es ) 2 = (Sx ™ + S 2 ) = (S x + S 2 ) 


dUpV 


, 2 /V 


Ui * U 
+ (Si* - S 2 2 )2 



u/ 


/ .2 / V A 2 

— + (Si - S a ) (— y 


4 (ei - ^ ^ - s 2 sgy - ( Sl - s 2 )= (Hi; 


dU 2 \ 2 _ >2 AiA 2 


t 


+ (Si 2 - S 2 2 )2 U V 


_ t(Sl + s 2 ) 2 (ll) 2 


(6a) 

(6b) 


In terms of equation (3), these become 


2 Si‘ 




U' V 

+ 2 

VU / U 2 



•\ 2 1 


u 


( 7 ) 


D 2 7 2 


S 2 2 r/ u «\ 


T LW 


- 2 


u c 


(-T1 

\u/ J 



to 2 7i +2 = (Sl + S 2 ) 2 (H_ j + ( Sl 2 _ S 2 2 )2 ^ 4 (Si - S 2 ) 


u 



2 ( VJ 


(8) 


( 9 ) 


to 2 7l -. 2 = (Sl - s 2 ) 2 (^y + (s x 2 - s 2 2 )2 ^ + (s x + s 2 ) 2 (zL 
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From equations ( 7 ) and ( 8 ), 


4r 5 


rA _ , fv 


s/ s/J 


U\ u/ V u 


and 


(n) 


Jf_ (7l + 7l) . §4^! 


fe)\(zi?l +2 (§X^)»4 (12) 

Xu/ \ u / J \ SxS 2 J u 2 


From equation (10) 


4D< . S x 2 + S 2 2 ./u ’\ 2 /v‘\a-| /Si 2 - 3 2 2 N u*v* 

(7x- 2 ) = - (- I + [— } ! + 2 (— — I 

SiS 2 L\x 


s x s 2 


x>/ 


V'u J ' 


S X S 2 J U” 


Subtracting equation ( 13 ) from equation (3.2) leaves 


, A1Y.A1Y 

c A u / \U ^ 
(13) 


4i>4_ 

S-.Sc 


(7i + 7s - 7i- £ ) = 2 


M - 

Au/ \U/ J 



(14) 


Adding equations (ll) and (l4) yie3.ds 


4 .P 2 r§ £ 

S X S 2 !-S x 


51 

7i + — 7 2 + 7i + 7 2 - 7 i- 2 

5 2 



>U 


\ 2 


By subtracting equation (l4) from equation (ll). 


4P 2 

SiS 2 


S 2 . \ 1 

— 7i + “ 7 2 - (7i + 7 2 ~ 7 1—2' 
-Si S 2 J 


-«*T 


and when equation (8) is subtracted from equation (7)> 

n 2 f 7i 72 V u,yt 

p / - p 


x 2 Sx 2 / u 2 


Lot (&f a o and + 

W 7i S 2 V 7i 


7 1—2 

7i y 


= a; then these three 


equations become 
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(15) 


(16) 


(17) 


And the correlation coefficient is 


1 - s 1 -- c r 1 ~ (l8) 

/u^/P 5 ^II TTT"a) (T"+l' ~a7 1 + V l“(rfrJ 

The relation 7 1+2 + 7i_ 2 = 2(7! + 72) gives alternate forms of the 
equations . 

In the case in which a is not 45°, 


(aif and 

\ uy / u >v v ' 2 

remain the same, hut 

1 — ) = (1 + £ - a) tan 2 a (19) 

\ u/ Si^ 


— — = D 2 (l - C) tan a (20) 

u 2 Si 2 

Equations (15), (19 ) , and (20) were not used in calculating the results 
presented in this paper. It was found that the conventional direct 
methods of measuring u’ and v ! gave better results. Equation (18) 
has the advantage of containing only ratios of the wire sensitivities 
and thermogalvanometer readings; thus it minimizes the possible errors 
involved in the readings. 


(£)* * * ST* (1 + 5 + a) 




= D 2 (1 + i - a) 


\ — j = d r 

Vu/ Si 2 


]i!zl = D 2ii d-n 

U 2 D Si 2 U * ' 
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APPENDIX B 

CALCULATION OF THE TURBULENT-SHEAR 
DISTRIBUTION FROM VELOCITY PROFILE 


Consider the "boundary— layer form of the Navier— Stokes equation 

u r + v r = i r (1) 

ox ay p oy 


Since it is known that the widths of the free turbulent mixing zones in- 
crease linearly with the axial distance, the following new independent 
variable can be introduced. 



If a stream function of the form \|r 
velocity components become 


' xFU) 


is chosen the two 


u = AF'(|) 

v = £ [|F*(i) - F(s) ] 


By substituting these forms in equation (l). 


-A 2 ^ F'F" + A 2 F" - (|F» - F) = ^ I G ' (0 

X X P x 


where r = BG(|) 


or 


-A 2 F"F = - BG* 
P 


G = - ^ / F"F di 

a B J 

r ^ 

f n 

a2- J 


+00 
'F d| 


A p +00 
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In order to be able to carry out the integration graphically, the 
form of the integral has to be changed. 


Assume F ! = 1 — f . Then 


I 


FF" d£ = 


-i-00 


FF® 


p 2 

- J f* as 


r £ +0 ° -,1 

« r (1 - t)(f a - f)ds + c ) - / (1 ~ f) 2 ds 


r. 


L 


,1 


= ' (1 — f)/ (1 - f)d| + c(l - f) - 


pi ni 

= -cf + j (f - f 2 )d| - f / (1 - f)d| 

+ CO 


f dl +f 2f 41 -/ f 2 dl 

-foo 


where c is the constant of integration. 
Hence, 


nrr 

pA 2 


r 

= -cf + / 




(f - f 2 )d| - f / (1 - f )d£ 


J 


From the boundary condition: 

t = 0, f = 1 for I = 

Therefore, 

7 —'00 —00 

c = j (f - f 2 )d! - f (1 - f)dt 

+00 

Thus, the dimensionless form of the turbulent shearing stress becomes 


- — • = -f / (f - f 2 )dt + f / (1 - f)dl + / (f - f 2 )dl 


PA 


+00 


+00 
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Figure 1.- Flow diagram of tlae two-dimensional tunnel. 
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Figs. 2,3 
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Figs. 4,5,6 




Figure 5.- Correlation figure at x * 54.2 
centimeters, y * 4.62. 



Figure 6.- Correlation figure at x * 54.3 
centimeters, y/$ * 3.72. 
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Figs. 7,8,9 
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Figure 7.- Correlation figure at x = 54.3 
centimeters, y/4 * 2.66. 



Figure 8.- Correlation figure at x = 54.3 
centimeters, y/4 * 2.05. 



Figure 9.- Correlation figure at x = 54.3 
centimeters, y/4 = 0.32. 
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rig«. 10,11 



Figure 10.- Correlation figure at x * 54.3 
centimeters, y/$ * 0.64, 



Figure 11.- Correlation figure at x * 54.3 
centimeters, y/$ ■ 2.86. 
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Fig. 13 MAGA TN No. 1257 
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figure 14.- Mean-velocity distributions (measured as compared with theoretical); d *= 12.0. 
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Figs. 16,18 
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Figure 16.- Lateral distribution of mean and fluctuating velocities and correlation 
coefficient; x * 30 centimeters. 
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Figure 18.- Lateral distribution of mean and fluctuating velocities and correlation 
coefficient; x » 75 centimeters. 
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Figure 17.— Lateral distribution of mean and fluctuating velocities and correlation 
coefficient; z ■ 54.3 centimeters. 


NACA TN No. 1257 


Fig®. 19 ,22 
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Figs. 20,21 



Figure 20 . — Long and short time exposure of 
correlation figure in the free- 
stream side of the let. 



Figure 21. - Long and short time exposure of 

correlation figure behind a screen. 


NATIONAL A DV I SON r 

COMMITTEE FOR AERONAUTICS 







NACA TN No. 1357 


Figs. 33,34 
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Figure 23.- Variation of exchange coefficient and mixing length across mixing rone. 
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Figure 24.- Ry distributions at different lateral positions; x ■ 30 centimeters. 
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Figs. 26,27 
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Figure 27.- R y distributions at different lateral positions; x ■ 54. 


centimeters. 


Figs. 28,29 
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Figure 28.- R y distributions at different lateral positions; x « 75 centimeters. 
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Figure 29.- R y distribution at x * 90 centimeters. 
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Fige. 30,33 
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Figure 33.- Comparison of R y and R 2 measurements; x =■ 75 centimeters, y/A - 0.61* 
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Figs. 31,32 



Figure 31.- 
y/* * 1.99 


* Long and short time exposure of R y 
correlation at x = 75 centimeters, 
with d = 0.027 centimeter. 




Figure 32 . 
y/$ = 1.99 


> Long and short time exposure of R y 
correlation at x = 75 centimeters, 
with d = 0.192 centimeter. 
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Figs. 34,35 
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Figure 34.- Ry dietribution at x « 30 centimeters, y/* - 0.47 indicating method of 
estimating the micro scale and scale of turbulence. 
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Figure 35.- Distribution of \ and L along the mixing tone. 
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Figs . 36,37 
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Figure 36.- Production, diffusion, and dissipation of turbulent energy across the 
mixing region. 
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